Thermal and electrical conductivity of Fermi pocket models of underdoped cuprates 
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Several models of the electronic spectrum in the pseudogap state of underdoped cuprates have been proposed 
to explain ARPES and STM measurements, which reveal only truncated Fermi pockets instead of a full metallic 
Fermi surface. We consider the transport properties expected of four physically distinct models, and calculate 
the thermal and electrical conductivity of the electronic quasiparticles. By proposing transport currents that 
reflect the close correspondence between quasiparticles on the Fermi pockets in the pseudogap and those near 
nodes in the superconducting state, we show that measurable transport coefficients provide stringent tests of 
pseudogap models. 

PACS numbers: 71.27.+a,72.10.Di,71.10.Ay 
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Angle-resolved photoemission (ARPES) and scanning tun- 
neling microscopy (STM) measurements on underdoped 
cuprates outside the superconducting dome have elucidated 
the electronic structure of the pseudogap statellHal- The large 
Fermi surface of overdoped systems is truncated to leave small 
arcs near nodal points. Beyond the arcs, quasiparticles are 
gapped, with the gap growing larger towards the antinodal re- 
gions. When the temperature is decreased below the super- 
conducting transition, the arcs themselves become gapped by 
a d- wave superconducting order parameterJ^-H] . 

Transport measurements on underdoped systems comple- 
ment the picture from ARPES and STM. The thermal k(T) and 
electrical cr(T) conductivity in the pseu dogap state transitions 
smoothly into the superconducting state 11914 1211. where the Bo- 
goliubov quasiparticles are well-understood. Wiedemann- 
Franz violations in the impurity- scattering regime, with the 
thermal conductivity exceeding the value expected from the 
Lorenz ratio k/(ctT) = Lq = (&g7r) 2 /3e 2 , are seen. Such viola- 
tions are familiar in the d- wave superconducting state where 
fermionic quasiparticles in the superconductor carry heat and 
charge at different velocities II13U 1511 . 

Several models of the pseudogap state, reviewed by Nor- 
man et al.[16], introduce coherent parts of the single electron 
Green's function G(k, co) and can account for features of the 
ARPES data. (The origin and properties of these models, as 
well as their respective success in capturing the ARPES and 
tunneling spectra, are discussed in detail in Ref. and refer- 
ences therein.) Each model Green's function describes a metal 
with two coherent quasiparticle bands that intersect near nodal 
points with an electron distributed between the two bands ac- 
cording to fractional k-dependent weights, just like in J- wave 
BCS theory. Slave boson treatments of Hubbard-models, as 
well as simple mean field approaches, have yielded Green's 
functions of this formfUM EH-El. 

The observable differences between the Green's functions 
for the candidate models are subtle, so the interpretation of 
additional experimental probes is crucial. In this article, we 
study the quasiparticle transport properties of representative 
models. Previous theoretical studies have obtained contra- 
dictory expressions for the quasiparticle electrical current- 
there is no obvious way to obtain a quasiparticle electrical 
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TABLE I: Four models of the pseudogap state. Each model has a 
single electron Green's function of the form Eq. [T]and accounts for 
some aspects of the ARPES spectra. They differ in the form of e k . 
The quantity e k is a tight-binding band energy expression that in- 
cludes several hopping terms while e 0k is restricted to nearest neigh- 
bour hopping. The abbreviations are defined in the text. 



current since its charge is not well-defined! 25- 29]. We pro- 
pose electrical and thermal currents that reflect the close cor- 
respondence between the observed quasiparticle transport in 
the pseudogap and J- wave superconducting states. Using this 
approach, we find violation of the Wiedemann-Franz law and 
differences in the transport properties of models that would be 
difficult to distinguish using ARPES spectra alone. 

Most of the single-electron Green's functions for the pseu- 
dogap models, reviewed in Ref. 16, can be written: 
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are the spectral weights of the '+' and '-' quasiparticle bands, 
E k = + A k their respective energies, and //k = (1/2) [6k + 
6 kl' ^ k ~ (1/2) [6k - 6 k ]. The pseudogap takes a simple J- wave 
form: Ak = Aol cos k x -cos k y \ while 6k and 6 k are tight-binding 
band energy expressions where the distinct form chosen for 6 k 
characterizes each model (see Table 1). 

To calculate the appropriate transport coefficients we start 
by assuming that we can write down a free-fermion Hamilto- 
nian with energies given by the poles of the Green's function: 
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where 7k,o-,± are quasiparticle operators. The corresponding 
(free quasiparticle) matrix Green's function G qp (k, co) is diag- 
onal, with respective elements: (co - E k+ )~ l and (co - E k _)~ l . 

For the model Green's functions above, electrons are frac- 
tionally distributed between bands. So quasiparticles will not 
carry charge at their band velocity and the scattering of elec- 
trons will couple the two quasiparticle bands (as is the case for 
superconductors To obtain simply heat and charge 

currents we transform to a basis of electrons, guided by the 
analogy with the J- wave superconductor. 

We thus assume a Bogoliubov transformation: 



U 



-w K . 



(4) 



which, when applied to the free-quasiparticle Green's function 
G qp (k, co) yields Eq. [T]as the [1 1] matrix component. Heat and 
charge currents are defined as for the quasiparticles in BCS, 
and transformed back into the quasiparticle basis. The thermal 
velocity is taken to be 
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where, in addition to the band velocities v+ and v_ there ap- 
pears a quantity v = -(^/E^dA^/dk + (A^/E^d^/dk. The 
latter mixes the contributions (to k(T)) of the two quasipar- 
ticle bands in the presence of a phenomenological scattering 
rate y. 

For the electric charge velocity we use: 
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where V/ = de^/dk and = de^/dk. 

We insert these expressions into the bubble diagrams for the 
thermal k/T and (quasiparticle) electrical conductivity <x: 
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where G R (k,co) includes self-energy effects associated 
with scattering (here we assume scattering is momentum- 
independent so vertex corrections can be ignored). 

Previous studies of the d-density-wave-, and BCS-derived 
models have disagreed J25I - E9I1 on the definition of the elec- 
trical current (Eq. [6]), particularly on whether to include off- 
diagonal elements associated with the velocity V2 = dA^/dk. 
If these elements are included, then y t h and \ e i are equal 
so Wiedemann-Franz violation does not occur ll27ll . The ob- 
servation of such violation and, more generally, the smooth 
evolution of quasiparticle transport across the superconduct- 
ing transition 111 Oil , supports the assumption that, like in BCS, 



the electric current does not include the off-diagonal terms. 
(Note that Eqs. [6] and [5] generalize the results for J- wave 
superconductors 113 On to cases where 6k ^ _6 k ) We now con- 
trast the transport properties of the models listed in Table 1 . 

The CDW model uses = 6k' where k' = k + Q and Q = 
(7r,7r). Ak couples electrons at different k-points in the full 
Brillouin zone, a quasiparticles is a combination of electrons 
at k and k + Q. The quasiparticle creation operator is: y^ ± = 
W+Ck ± W-Ck+q and the appropriate electric velocities are then 
V/ = de^/dk and = de^ /dk\ which for nearest-neighbour 
hopping gives = -\f. The two electronic components of 
the quasiparticle carry current in opposite directions. 

The BCS Green's function is obtained by taking = -6k. 
A quasiparticle is built from the coupling of an electron at mo- 
mentum k and a hole of momentum -k and the quasiparticle 
are y^+ = W+c^± W-c[_ k . Since both the velocity and charge 
of the quasiparticle 's constituents differ in sign, they carry 
equal current so electric transport velocities satisfy v/ = v^. 
The sign change relative to the CDW case has crucial conse- 
quences for the Lorenz ratio, as seen below. (Henceforth, we 
refer to the BCS-like relationship between the electron and 
hole velocities: v/ = as electron-hole coupling, and to 
the CDW-like relationship = -v/ as electro-electron cou- 
pling.) 

For nearest neighbour hopping, the CDW and BCS models 
can be obtained by appropriate choice of parameters in an- 
other model discussed in Ref. LUi the energy-displaced node 
model (EDN). The band energies satisfy = -6k and the pa- 
rameter yUk = //o is constant. A finite jiq provides for a finite 
zero-energy density of states (i.e. a Fermi pocket with a radius 
proportional to //o) as illustrated in the inset of Fig 1. At low 
temperatures, ksT « max(yUo,y) « Aq we find: 
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and, taking = -V/ as appropriate for the CDW model, we 
have 



cr = o- (k/ko) 



(10) 



where Kq/T = L [l + a 2 ] 2 , where a = v^/vf and <x = 
(e/7t) 2 a~ l . A similar result is obtained in Ref. |2C|. In the limit 
y » yuo we recover the universal values for transport coeffi- 
cients. In the limit yUo » y we find (k/T) — > (Ko/T)[njuo/4y] 
and cr — > <xo0uo7r/4y). So the conductivities scale with the ra- 
dius of the Fermi pocket yUo and with the quasiparticle lifetime 
1 ly as expected for a simple metal. Note that (L/Lo) = [1 +a 2 ] 
is independent of jiq /y, so there is a constant enhancement of 
the heat conductivity relative to the Wiedemann-Franz predic- 
tion. In Fig. 1 we show the conductivity in the EDN model 
with nearest neighbour hopping as a function of y for several 
values of yUo. Evident is the crossover from behaviour at small 
yuo to linear y~ ^dependence. 

While none of the models discussed in Ref. [l6| consid- 
ered a finite j^o in the presence of an electron-hole coupling, 
it is simple and interesting to extend our analysis to such 
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FIG. 1 : Fermi pockets, nodes and universal transport. Several model 
Green's functions of the pseudogap regime have quasiparticle dis- 
persions resembling those in d- wave superconductors. Near a node 
(indicated by the circle appearing along the zone diagonal on the 
normal state Fermi surface shown in the upper left) there are two in- 
tersecting quasiparticle bands, as illustrated. A key feature of some 
of the models (including the EDN model, plotted here) is an energy 
shift, /io, of the Fermi level from the node, which gives rise to a Fermi 
pocket. The curves show the electrical conductivity cr in units of the 
universal value cr plotted versus scattering rate y/A for different 
values of yu . When A /y is large and quasiparticles are restricted to 
nodes, the value of /i determines whether quasiparticle transport be- 
haves as in a normal metal (with conductivity proportional to y _1 ) or 
as in a d- wave superconductor (with universal conductivity cr = cr ). 
For small y/A , quasiparticles are excited everywhere on the Fermi 
surface and the parameter jd Q plays no role. 



models. Using \f = \' the thermal conductivity k, Eq. [9l 
is unchanged but the electrical conductivity becomes cr = 



'(f). 



which again gives cr in the jdo/y « 1 

universal limit. However, in the fio/7 >> 1 limit, gives a value 
twice as large as the corresponding limit of Eq. [TUJ (An anal- 
ogous factor of two was found |28|] at finite temperature when 
fio = 0.) So, for electron-hole couplings, L/Lq is strongly 
y-dependent and can be either smaller or larger than one. 

We now contrast the ARPES and transport predictions of 
the models discussed above. In the CDW model, the pre- 
dicted ARPES peaks trace out the elliptical Fermi pocket with 
axes of length jio/vf and aj^o/vf. (Because of the spectral 
weight, a portion of the back side of the pocket would be ob- 



scured, leaving only an arc.) The corresponding transport pre- 
diction is that the conductivity should scale with [jiq /a](l +a 2 ) 
whereas the Lorenz ratio should be L = Lo(l + a 2 ). If jiq and 
v//v2 changed differently with doping then one could impose 
a strong consistency check on the doping-dependence of the 
ARPES spectrum. 

In the BCS model there are no pockets but there do appear 
arcs in the presence of a scattering rate y since the poles at en- 
ergies ±£k are smeared together giving an arc of zero-energy 
spectral peaks along 6k = 0. The length of the arc is pro- 
portional to y/v2 but the transport coefficents are universal, so 
changes in the scattering rate for a given size of the pseudo- 
gap would change the length of the arc without affecting the 
conductivity. 

As a concrete example of a possible doping evolution of 
the conductivities we consider the YRZ model B2 1I423H ob- 
tained by taking = -6ok = -2t[cosk x + cosk y ], and 
6k = -2t[cosk x + cosk y ] - 4t' [cos k x cos k y ] - 2£"[cos2& x + 
cos 2k y ] + yuo where all the hopping coefficients and jiq are de- 
pendent on doping in a precise manner prescribed by YRZ 
(f,f and yuo vanish at half-filling and increase linearly in 
magnitude with doping away from half-filling). The result- 
ing doping-dependence of the Fermi pockets is illustrated in 
the sketches on Figs. 2 and 3. In addition, the weight of 
the quasiparticle is zero at half-filling and increases approx- 
imately linearly. 

The electric current velocities are given by v/ = Vk6k and 
y'j. = +Vk6 k where the sign corresponds to the electron- 
electron (i.e. akin to CDW) and electron-hole (i.e. like BCS) 
coupling possibilities, respectively. 

If Fig. 2 we show the result for electron-hole coupling with 
the conductivity plotted in the main panel and the Lorenz ra- 
tio shown in the inset. The rapid increase of the conductivity 
with doping results from the growth of the Fermi pocket and 
the quasiparticle weight. The Lorenz ratio is strongly dop- 
ing dependent, especially when the scattering rate is low. For 
small y, L/Lq attains large values at low doping (where the 
magnitude of the pseudogap is large), dips below one at in- 
termediate doping and finally approaches one at large doping 
as the pseudogap vanishes. In Figure 3, we contrast the case 
of electron-hole and electron-electron coupling, again using 
the doping dependence prescribed by the YRZ model. The 
biggest difference between the two pairing scenarios is the 
strong y-dependence seen in the electron-hole coupling but 
not the electron-electron coupling. 

The existing thermal and electrical transport data on non- 
superconducting, underdoped cuprates supports the ARPES 
though may not yet rule out any models discussedll9l-[l2ll. The 
quasiparticle contribution to conduction decreases with un- 
derdoping and large Lorenz values (L/Lq > 3) have been 
seen in the field induced normal state in metallic, underdoped 
samples II 1211 . The Lorenz ratio is larger for dirtier samples 
(where the ratio of low and high temperature resistivies was 
used to estimate sample purity). These results are in qualita- 
tive agreement with Fig. 3 and the dependence of the Lorenz 
ratio on scattering rate suggest that quasiparticles are formed 
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FIG. 2: Doping-dependence of electrical and thermal conductivi- 
ties. The electrical conductivity <t/o- is plotted versus doping 6 
in the main panel, assuming the YRZ-evolution of electronic band 
parameters 1 22, 23], and quasiparticle weight, for several values of 
the scattering rate y (in units of bare, doping-independent nearest- 
neighbour hopping parameter t). The change of the Fermi pockets are 
illustrated by the small figures. In the inset, the Lorenz ratio L/L , 
where L is the Lorenz number, is plotted versus doping (showing 
deviations from the Wiedemann-Franz law L = L ). Values of L/L 
significantly greater than one can occur at small 6, owing to a large 
contribution to heat transport associated with momentum-variation 
of the pseudogap. In clean systems (i.e. for small y) with an electron- 
hole Bogoliubov transformation assumed (see text), Lorenz values of 
less than one are also possible. The qualitative doping dependence 
of the Lorenz ratio depends on model parameters, and can be used 
to distinguish candidate models (especially when the individual con- 
ductivities are not easy to interpret). 



from electron-hole coupling. However insulating behaviour 
has been reported that, in some studies, effects charge but not 
heat transport! 9], which would indicate a more radical depar- 
ture between heat and charge currents than proposed above. 

In conclusion, transport coefficients of the model Green's 
functions, which describe coherent quasiparticles near nodes, 
offer tests of the models beyond those provided by the ARPES 
spectra that motivated them. To combine the benefits of both 
probes, we have provided a plausible conjecture for the heat 
and charge currents that captures the observed connection be- 
tween the pseudogap and the d-wave superconducting quasi- 
particles. 

This work was supported by Australian Research Council 
Discovery Projects DP1094395 and DP0710617. We thank 



FIG. 3: Scattering dependence of the Lorenz ratio in Fermi pocket 
models: contrasting electron-hole and electron-electron coupling. 
The Lorenz ratio L/L is plotted versus the scattering rate y for sev- 
eral values of doping 6. The doping dependence of all parameters is 
the same as in Fig. 2. If the quasiparticles are formed from a combi- 
nation of an electron and a hole (resulting in dashed curves in plot) 
then the Lorenz ratio is strongly dependent on scattering rate. If they 
are formed from a pair of electrons at different momenta (resulting 
in solid curves) then the Lorenz ratio is independent of scattering 
rate. The qualitative behaviour is significantly different between the 
electron-hole and electron-electron pairing scenarios (whereas both 
can give rise to the same single-electron spectral function as seen in 
ARPES data). Experimental transport data on far-underdoped sys- 
tems suggest a strong dependence of the Lorenz ratio on scattering 
rate in support of the electron-hole picture. 
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